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The heat capacity of the binary liquid mixture triethylamine-water has been measured near its lower 
critical consolute point using a scanning, adiabatic calorimeter. Two data runs are analyzed to 
provide heat capacity and enthalpy data that are fitted by equations with background terms and a 
critical term that includes correction to scaling. The critical exponent a was determined to be 
0. 107±0.006, consistent with theoretical predictions. When a was fixed at 0. 1 1 to determine various 
amplitudes consistently, our values of A + and A~ agreed with a previous heat capacity 
measurement, but the value of A + was inconsistent with values determined by density or refractive 
index measurements. While our value for the amplitude ratio A + M“=0.56±0.02 was consistent 
with other recent experimental determinations in binary liquid mixtures, it was slightly larger than 
either theoretical predictions or recent experimental values in liquid-vapor systems. The correction 
to scaling amplitude ratio D + /D~= 0.5±0.1 was half of that predicted. As a result of several more 
precise theoretical calculations and experimental determinations, the two-scale-factor universality 
ratio X, which we found to be 0.019±0.003, now is consistent among experiments and theories. A 
new universal amplitude ratio R H involving the amplitudes for the specific heat was tested. Our 
determination of R B ^ = -0.5 ± 0.1 and R Ba ~ ~ 11 ± 0.1 is smallerin magnitude than predicted and 

is the first such determination in a binary fluid mixture. © 1996 American Institute of Physics 
[S002 1-9606(96)50820-2] ' 


INTRODUCTION 

Critical point phenomena have been studied over a long 
period of time with substantial recent effort toward modeling 
complex systems both near a critical point and over an ex- 
tended region in thermodynamic space. Large density or con- 
centration fluctuations near a system’s critical point effec- 
tively mask the identity of the system and produce universal 
phenomena which have been well studied in simple liquid- 
vapor and liquid— liquid systems. 1-3 Such systems have pro- 
vided useful model systems to test theoretical predictions 
which can then be extended to more complicated systems. 
Along various thermodynamic paths, several quantities ex- 
hibit a simple power-law dependence close to the critical 
point. The critical exponents describing these relationships 
are universal and should depend only on a universality class 
determined by the order-parameter and spatial dimensional- 
ity of the system. Binary fluid mixtures, liquid— vapor sys- 
tems, and uniaxial ferromagnetic materials are thought 1 ' 1 to 
belong to the three dimensional, Ising model which, in a 
renormalization group (RG) context, predicts 4 values. for the 
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critical exponents. These values, for the most part, agree very 
well with experimental results. 1-3 Several universal ampli- 
tude ratios have also been predicted and many have also 
been experimentally tested with mixed results. 5 A particu- 
larly powerful observation is that only two critical exponents 
are linearly independent, and that the leading amplitudes are 
interrelated using only two scale factors. 1-3 Thus the univer- 
sality of the exponents could be used with two experiments 
to determine all the leading critical behavior of a given sys- 
tem. 

The heat capacity provides a delicate probe of the sys- 
tem near a critical point and can determine essential ampli- 
tude and exponent values. In particular, a precise measure- 
ment of the heat capacity of the binary fluid mixture 
triethylamine and water will help settle a dispute in the lit- 
erature about the amplitude in the one-phase region and its 
effect on the universal ratio X, measure the universal ratio of 
the amplitudes above and below the critical temperature, and 
give the first attempt in binary liquid mixtures at using a new 
“universal” amplitude relation R^ to reduce the number of 
fitting parameters while testing theoretical predictions. 

The heat capacity has a weak divergence near the critical 
point that is governed by the critical exponent a, has a criti- 
cal contribution B c to the background heat capacity, and has 
correction to scaling terms that extend the theoretical de- 
scription further from the critical point. 1 ’ 6 The critical point 
influences the system properties over a global region; 7 8 how- 
ever, this experiment measures behavior near the critical 
point and we will use the simpler formalism that is appropri- 
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ate there. While a vapor-liquid system measured along the 
critical isochore would have a weak divergence in the heat 


capacity at constant volume, a binary fluid mixture near its 
critical consolute point will have the same type of divergence 
in its heat capacity at constant pressure when measured along 
a path of constant, critical composition x~x c ; 3 

A + 

C P x=(Bt> g + B c ) + Et+ — \t\- a (l+D + \t\*' 

-» ) one phase, (j a ) 

A - 

C P x =z (B hg + B c ) + Et+ — M~ a (l+D'|r| A i 

H ) two phase, (lb) 

where C px is the temperature dependent heat capacity at con- 
stant pressure and at the critical composition, B hg is the 
background heat capacity far from the critical temperature T c 
in the one-phase region, t — (T—T c )fT c is the reduced tem- 
perature, A is the amplitude of the leading divergence, 
D + is the amplitude of the first correction to scaling term, 
and A! is the universal exponent for the correction term. The 
one and two phase amplitudes are denoted by the -I- and - 
superscript respectively; the background terms and the criti- 
cal exponents are predicted 6 to be the same above and below 
the critical point. A linear background term Et arises from 
the regular part of the free energy and hence should be the 
same above and below the critical point. The critical expo- 
nents a and A, are predicted 9 to be 0.110±0.003 and 0.51 
±0.03, respectively. 

The critical part of the background heat capacity has an 
interesting connection with some of the other amplitudes 
governing the heat capacity. Bagnuls and Bervillier have 
predicted 10 a universal amplitude ratio defined by 

D Cr J 

„ i4 ± |D ± |“ M ' 

^B c * ® 

with predicted values of R ^ = -0.708 15 for the ratio in- 
volving one-phase" amplitudes A and D, and R B 

= - 1.334± 0.044 for the two-phase 6 region. One of their 
motivations was to allow an experimental test of the theory 
using fewer adjustable parameters. Others 7812 have noted 
that Rg^ is not universal in the same sense as the universality 
of the critical exponents. Bagnuls and Bervillier 13 have 
agreed, but noted that Rg may be universal for a subclass of 
systems, and that the amplitude ratio D* ID (predicted 6 to 
be 0.96±0.25) is similarly not universal. 

Other amplitude ratios are predicted to be true universal 
quantities that should be the same for all systems in a given 
universality class (three-dimensional Ising for bulk, binary 
liquid mixtures). The ratio of the leading, singular, heat ca- 
pacity amplitude in the one-phase region to that in the two- 
phase region is now predicted to be A + /A ~ =0,541 ±0.014 
by RG 6 and 0.523±0.009 by series 9 calculation. One can 
also relate A + to the amplitude §, of the correlation length in 
the one-phase region defined by £=&/-*' using two-scale- 


factor universality: 9 X = A + £/k B , where k B is Boltzmann’s 
constant, and v is another universal critical exponent whose 
value 4 is about 0.63. The value of X is predicted from RG 14 
to be X=0.019 66±0.000 17 and from series 9 to be 0 0188 
±0.000 15. 

These predictions have been tested in other binary fluid 
mixtures with mixed results. 1 3 Recent experiments measur- 
ing the heat capacity of binary fluid mixtures have given 
A I A in the range 0.56 to 0.71; 13 18 a value consistently 
larger than the 0.52 value determined in liquid-vapor 
systems. 19 To determine the universal quantity X , the value 
of the correlation length amplitude must be determined in 
a separate experiment. For the system triethylamine-water, 
the reported values of & have varied from 0. 10 nm (Ref. 15) 
to 0.125 nm (Ref. 20), with the most recent determination 
being 0.1 10±0.005 nm. 21 Since £ 0 enters as the cube in X, 
the error in usually dominates the error in X. A number of 
values of X, ranging from 0.017 to 0.024, have been deter- 
mined using a variety of values for £ 0 and A + just for the 
system triethylamine and water. Recent determinations in 
other binary fluid systems have found X to be 0.019±0.004 
in triethylamine and heavy water, 17 and 0.020±0.002 in 
3-methylpentane and nitroethane. 18 

The value of A + for the triethylamine-water system re- 
ported here has had a number of determinations using vari- 
ous techniques. Thoen et a/. 15 measured the heat capacity 
using a step technique in an adiabatic calorimeter which did 
not stir the sample. They found A + to be 0.23±0.05 
J/(cm 3 K). Three other publications 17-21,22 have reanalyzed 
those data and determined A + to be between 0.23 and 0.36 
J/(cm K) (see Table I). One can also determine the ampli- 
tude A + from precise density measurements, along with the 
slope of the critical temperature dependence on pressure. 
Such determinations 5-2 1-22 have given A + between 0.1 19 and 

0. 35 J/(cm 3 K) as shown in Table I. Another possibility is to 
measure the refractive index and then calculate the density; 
however, the refractive index itself may have an intrinsic 
anomaly which is an issue not yet resolved. 23-24 Values of A + 
determined by refractive index tend to be small, while values 
determined by density tend to be large, as is evident in Table 

1. One of the motivations of this experiment was to resolve 
this controversy by direct measurement. 


EXPERIMENT 

The heat capacity of a critical mixture of triethylamine 
and water was measured using a scanning adiabatic calorim- 
eter similar to that used by Thoen et al . 15-17 for studies in this 
same system and in triethylamine and heavy water. Our calo- 
rimeter improves upon Thoen’s 15 earlier experiment on this 
system by including stirring of the sample, operating in a 
scanning mode, and using computer data acquisition and 
control. An important advantage over ac calorimetry 18 is our 
ability to measure the background heat capacity, which is 
important in theoretical constructions. A very slow tempera- 
ture scan through the critical point is essential to allow the 
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TABLE I. Amplitude of the specific heat anomaly in the one (A + ) and two 
M“) Phase region for triethylamme-water (TEA-H 2 0) as found by differ- 
ent techniques. The refractive index n and density p anomalies allow A* to 
be calculated using thermodynamic relations. The amplitudes determined by 
analyzing prior heat capacity data (Ref. 15) depend on the choice of critical 
temperature and parameterization. A similar system, triethylamine-heavy 
water (TEA-D 2 0), is provided for comparison. All errors are one standard 
deviation estimates as given by the authors. 



A + (J/cm’ K) 

A ~(J/cm 1 K) 

System/approach: 

tea-h 2 o 

n 

f 0. 1 16± 0.005 a 
| 0. 15±0.01 b 
[0.271±0.007 C 


P 

r 0.1 19± 0.009 d 
0.35±0.02 c 
. 0.267±0.002 c 


C r 

r0.23±0.05 f 

0.23±0.02 g 

0.246±0.005 h 

0.36±0.02' 

0.201 ± 0.006 k 

0.33±0.06 f 

0.358 ±0.005 k 

TEA-D.O 

c P 

0.215±0.01 J 

0.378 ±0.004' 

“Reference 20. 

‘’Reference 23. 

Reference 21. 
d Reference 5. 

Reference 22. 

Reference 15. 

8 ReanaIysis of Ref. 15 in Ref. 17. 
h Reanalysis of Ref. 15 in Ref. 21. 
'Reanalysis of Ref. 15 in Ref. 22. 


^Reference 17. 
k Our results in bold. 


system to stay in quasistatic equilibrium throughout the ex- 
periment and thus considerable time and effort goes into the 
design and construction of the thermostat. 

A critical composition sample of triethylamine and water 
was used. The triethylamine was 99% pure from Kodak and 
was further distilled over KOH with the center fraction saved 
under a dry nitrogen atmosphere in a Teflon bottle. The wa- 
ter was from E-M Science: “Omni-Sol v glass distilled’ ’ 
which was frozen and vacuum pumped to remove dissolved 
gases before storage under nitrogen. The mixture used in this 
experiment was 31.7% by weight of triethylamine in water 
with a total mass of 16.015 g which is a volume of 17.20 mL 
at the critical point. This concentration of triethylamine and 
water is the most recently determined critical concentration 25 
and is close to the 32.1% used by other investigators of this 
system. 5rl5J7_,8 ’ 20 “ i ' 4 ' 26 ’ 27 The sample preparation was done 
under a nitrogen atmosphere and a small nitrogen and vapor 
bubble was left above the fluids to keep the pressure close to 
one atmosphere. This system exhibits a lower consolute 
point around 18.2 °C and hence is one-phase below the criti- 
cal point and two-phase above. Triethylamine is a strong 
base which is very reactive and requires special containment. 

The mixture was sealed in a cylindrical, gold plated, 



FIG. 1. Our adiabatic calorimeter is a set of temperature controlled, concen- 
tric cylindrical shells. The fluids F are sealed in a gold-plated cell and are 
stirred by a stirbar coupled to a stir motor SM attached to the outer stage O. 
A passive radiation stage RA surrounds the cell and is contained within the 
reference stage RE which is actively temperature controlled to be at the 
same temperature as the cell. Stages are separated by nylon standoffs N. The 
inner stage I and vacuum stage V are also temperature controlled. A vacion 
pump VP reduces heat conduction losses by maintaining the pressure inside 
the thermostat to 10 6 torr while electrical signals pass through a wire 
feedthrough WF. 

oxygen-free high conductivity copper cell with a Kalrez (Du- 
Pont perfiuoroelastomer) o ring. The fluids were stirred with 
a Teflon coated stir bar which was magnetically coupled to a 
small motor attached to the outer stage of the thermostat as 
described below. Stirring the mixture is essential to achieve 
thermal equilibrium in the two-phase region. 17,1821,25 The 
fluids did not react with the cell as evidenced by a critical 
temperature which varied less than 2 mK over a seven month 
period, which is in contrast to a 1 mK/day drift when these 
fluids are sealed in glass. 21,26 The cell is heated with a con- 
stant current (to 0.01%) flowing through manganin wire 
wrapped in grooves on the outside of the cell body. 

The cell is surrounded with a nested set of cylindrical, 
nickel-plated, copper stages forming a thermostat. A passive 
radiation shield surrounds the cell and is in turn encom- 
passed by a reference stage that is actively temperature con- 
trolled so as to be at the same temperature as the cell 
throughout the experiment. Two additional actively con- 
trolled stages are maintained approximately 0.1 K below the 
next inner stage. The outermost of these has the stir motor 
attached to it. All the stages are attached to each other with 
hollow nylon stand-offs. This assembly is placed in a 
vacuum stage that has an ambient temperature determined by 
a Lauda RT-7 bath that circulates water through copper tub- 
ing soldered to the outside walls and lid, and is evacuated to 
1 ^atorr using a vacion pump. The thermostat is depicted in 
Fig. 1. 


J. Chem. Phys., Vol. 104, No. 20, 22 May 1996 





Flewelling et at.: Heat capacity anomaly in (CH 3 CH 2 )3N-H 2 0 


8051 


The temperature of the stages is determined by ther- 
mistors embedded through the lids and into the side wall of 
each stage. The cell and reference stage use matched, “ultra- 
stable.” Thermometries thermistors while the other stages 
use Yellow Springs thermistors. All of the thermistors are 
wired in series with a 100 kfl Vishay standard resistor and a 
10 /zA constant current source. The resistance of each ther- 
mistor is determined by measuring the voltage across the 
standard resistor and across each thermistor using a 7 j digit 
Keithley 2001 multimeter. The power into the cell thermistor 
is about 2 /zW which is a negligible 0.05% of the power 
provided by the heater. This small power in the thermistor 
means that its self-heating is less than the temperature repro- 
ducibility of 0.2 mK. Each thermistor has been calibrated 
relative to a Guildline platinum resistance thermometer, 
which provides an accuracy of 10 mK in absolute tempera- 
ture and 0.3 mK in relative temperature. 

The digital multimeter measures the thermistor and stan- 
dard resistor voltages and the voltage applied to the cell 
heater and reports each to a Macintosh Ilci computer using a 
GPIB bus. The computer data acquisition and control soft- 
ware is written in the programming language Lab VIEW and 
uses a software PID controlling algorithm to maintain the 
reference stage temperature within 1 mK rms of the cell tem- 
perature. Two Kepco GPIB, digital to analog converters are 
used to set the heater voltages on the reference, inner, and 
outer stages. The program measures, controls and records 
temperatures with a cycle time of 12 s. The cell heater re- 
ceives 4 mW of power that causes the cell temperature to 
drift about 30 /zK/s so that the cell temperature is effectively 
constant during each data cycle. 

A subsequent calibration run using only water deter- 
mined a number of experimentally important properties of 
the calorimeter. It was found that the power delivered to the 
fluids by the stirrer was less than 0. 1% of the power supplied 
by the heater and thus negligible in this experiment. Using 
the known specific heat of water, we found that the empty 
cell contributes 20.4±0.1 J/K to the total heat capacity. 
Small heat leaks between the cell and other stages were sys- 
tematically quantified so corrections could be applied to the 
raw data on the critical mixture. 

Two data runs are reported here on a critical composi- 
tion of triethylamine and water between 14.5 °C and 19.5 °C. 
Each data run takes almost two full days and entails over ten 
thousand sets of stage temperatures, cell temperature, and 
cell heater voltage. Early data runs gave a rounded anomaly 
and very noisy two-phase data as others 1718 have reported 
when stirring was not effective. The principle raw data are 
the variation in cell temperature with time, which can be 
analyzed in two ways to determine the amplitudes and expo- 
nents associated with the heat capacity. The first method is to 
determine the rate of change of the cell temperature dT ce]l /dt 
over some range of temperatures by doing a straight line, 
least-squares fit and then calculate the fluid’s average heat 
capacity at the average temperature directly: 


V 2 

— -/r J (7' ce , ,-r f ) 


" (dT cen /dt)V Q 


[J/(cm 3 K)], 


(3) 


where K, is the thermal conductivity between the cell at tem- 
perature T cM and stage / at temperature 7, , V is the average, 
measured voltage across the cell heater of constant resistance 
/?, Vq is the fluid volume, and the empty cell heat capacity 
has been subtracted. The only significant heat loss was be- 
tween the cell and reference stage for which we can correct 
using our calibration data on water. The number of points 
used to determine the cell temperature’s rate of change needs 
to be large enough to minimize the effect of the temperature 
resolution yet not so large as to average over a region where 
the heat capacity changes significantly. We have chosen to 
use 200, 100, 50, and 20 points to determine the slope 
dT ce \\/dt using the larger number of points furthest from the 
critical point. The heat capacities so determined from the two 
runs are listed in Table II and illustrated in Fig. 2. The errors 
represent the random fluctuations in the raw data and the 
propagated error from determining the slope of the cell tem- 
perature. It should be noted that the values of the heat capac- 
ity determined in this way are at a saturated vapor pressure, 
but Thoen et al.' 5 has demonstrated that the correction to 
constant pressure and composition is negligible and thus our 
values of the heat capacity are just C px . 

An alternative determination of the heat capacity param- 
eters is to analyze the enthalpy as a function of the cell 
temperature. The enthalpy H of a data run is calculated by 
adding the energy to the cell from the heater and subtracting 
the heat loss to other stages over each time interval At. 




y-Wen-T,), 


At, /IV 


(4) 


The enthalpy is then related to the heat capacity by a simple 
integral over temperature T : 


H = H ° + lr C "' dT (5) 

The calculated enthalpy divided by the critical temperature 
T c is illustrated in Fig. 3 for run 2. Using the enthalpy has 
the advantage of analyzing the entire data set, but the disad- 
vantages of accumulated error (resulting in systematic off- 
sets) and a (I —a) cusp instead of a divergence at the critical 
point. We will analyze our data using both techniques in the 
next section. 


ANALYSIS AND INTERPRETATION 

A weighted, nonlinear, least squares routine was used to 
fit the equations to the data by finding the best set of param- 
eters that minimize" 8 the reduced chi square ^IN. Using our 
own functions in the commercial program Igor Pro, 29 we 
could allow some parameters to be free and force others to 
take on certain values as will be described below. This pro- 
gram allows each data point to be weighted according to the 
error in the dependent variable and determines errors on the 
parameters that are related to the diagonal elements of the 
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TABLE II. Heat capacity data for two runs at the critical composition in triethylamine and water. Temperatures 
are in °C and the critical temperature is 18.225±0.001 for run 1 and 18.227±0.001 for run 2 SC is the 
propagated error in C p . p 


Temperature 

run 1 


C, 

J/(cm 3 K) 


SC 

J/(cm 3 K) 


Temperature 


J/(cm 3 K) 


SC 

J/(cm 3 K) 


run 2 


15.0121 

4.08 

15.2161 

4.10 

15.4293 

4.13 

15.6448 

4.15 

15.8616 

4.18 

16.0773 

4.21 

16.2945 

4.23 

16.5134 

4.27 

16.7328 

4.32 

16.9527 

4.37 

17.1148 

4.42 

17.2107 

4.46 

17.3052 

4.48 

17.3989 

4.52 

17.4922 

4.57 

17.5844 

4.63 

17.6762 

4.69 

17.7443 

4.74 

17.7878 

4.78 

17.8309 

4.84 

17.874 

4.88 

17.9143 

4.94 

17.9533 

4.99 

17.9919 

5.06 

18.0301 

5.13 

18.0678 

5.24 

18.1048 

5.35 

18.1411 

5.51 

18.1763 

5.76 

18.185 

5.85 

18.2018 

6.05 

18.2177 

6.80 

18.2222 

7.32 

18.2318 

10.02 

18.2361 

9.56 

18.2402 

11.66 

18.2444 

10.03 

18.2524 

9.89 

18.2644 

9.90 

18.2766 

9.70 

18.2891 

9.25 

18.3021 

9.20 

18.3086 

9.03 

18.3353 

8.67 

18.3623 

8.64 

18.3895 

8.47 

18.418 

8.29 

18.4464 

8.28 

18.4754 

8.02 

18.5047 

8.01 

18.5333 

7.91 

18.5614 

7.83 

18.589 

7.74 

18.616 

7.72 

18.6431 

7.64 

18.7117 

7.51 

18.7673 

7.44 

18.8234 

7.35 

18.8799 

7.30 

18.9368 

7.24 

18.9942 

7.18 

19.0519 

7.12 


0.01 

15.1085 

0.01 

15.3136 

0.01 

15.5206 

0.01 

15.7263 

0.01 

15.9334 

0.01 

16.1414 

0.01 

16.3504 

0.01 

16.5514 

0.01 

16.729 

0.01 

16.904 

0.02 

17.0359 

0.02 

17.1222 

0.02 

17.2106 

0.02 

17.2984 

0.02 

17.3831 

0.02 

17.4695 

0.02 

17.5551 

0.03 

17.6389 

0.03 

17.7006 

0.03 

17.7434 

0.03 

17.7851 

0.03 

17.8264 

0.03 

17.8668 

0.03 

17.9069 

0.03 

17.9479 

0.03 

17.9872 

0.03 

18.0254 

0.03 

18.063 

0.03 

18.1 

0.08 

18.1363 

0.08 

18.173 

0.09 

18.182 

0.23 

18.2004 

0.30 

18.2177 

0.29 

18.2223 

0.35 

18.2319 

0.30 

18.2361 

0.12 

18.2402 

0.12 

18.2443 

0.12 

18.2484 

0.12 

18.2514 

0.12 

18.2619 

0.06 

18.2728 

0.06 

18.284 

0.06 

18.3014 

0.06 

18.3252 

0.06 

18.3495 

0.06 

18.374 

0.06 

18.3989 

0.06 

18.4243 

0.06 

18.4497 

0.06 

18.4757 

0.06 

18.5018 

0.05 

18.528 

0.05 

18.5545 

0.03 

18.5812 

0.03 

18.5947 

0.03 

18.6491 

0.03 

18.7042 

0.03 

18.7602 

0.03 

18.8169 

0.03 

18.8741 


4.02 

0.01 

4.02 

0.01 

4.06 

0.01 

4.07 

0.01 

4.09 

0.01 

4.12 

0.01 

4.16 

0.01 

4.21 

0.01 

4.25 

0.01 

4.30 

0.01 

4.34 

0.02 

4.37 

0.02 

4.40 

0.02 

4.44 

0.02 

4.48 

0.02 

4.57 

0.02 

4.60 

0.02 

4.63 

0.02 

4.68 

0.03 

4.71 

0.03 

4.76 

0.03 

4.79 

0.03 

4.86 

0.03 

4.92 

0.03 

4.97 

0.03 

5.03 

0.03 

5.12 

0.03 

5.21 

0.03 

5.34 

0.03 

5.50 

0.04 

5.77 

0.04 

5.88 

0.08 

6.05 

0.08 

6.58 

0.08 

7.04 

0.20 

10.53 

0.31 

10.14 

0.30 

11.36 

0.34 

10.95 

0.33 

10.59 

0.32 

10.51 

0.13 

9.98 

0.13 

9.62 

0.12 

9.32 

0.12 

9.01 

0.06 

8.82 

0.06 

8.61 

0.06 

8.58 

0.06 

8.36 

0.06 

8.25 

0.06 

8.10 

0.05 

8.06 

0.05 

8.02 

0.05 

7.95 

0.05 

7.83 

0.05 

7.78 

0.05 

7.70 

0.03 

7.58 

0.03 

7.48 

0.03 

7.35 

0.03 

7.25 

0.03 

7.19 

0.03 
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TABLE II. ( Continued .) 


c P 

sc 

Temperature 

J/(cm 3 K) 

J/(cm‘ K) 

19.1107 

7.10 

0.03 

19.1692 

7.07 

0.03 

19.2277 

7.02 

0.03 

19.2865 

6.99 

0.03 

19.3456 

6.96 

0.03 

19.4049 

6.90 

0.03 

19.4645 

6.89 

0.03 



C r 

SC 

Temperature 

J/(cm' K) 

J/(cm f K) 

18.9316 

7.16 

0.03 

18.9892 

7.13 

0.03 

19.0472 

7.08 

0.03 

19.1055 

7.03 

0.03 

19.1642 

6.97 

0.03 

19.2232 

6.94 

0.03 

19.2825 

6.92 

0.03 

19.3419 

6.89 

0.03 

19.4015 

6.85 

0.03 

19.4615 

6.82 

0.03 


typical error matrix. 28 While this is adequate for equations 
whose parameters enter linearly, errors of nonlinear param- 
eters should be determined by error ellipses. 28 In tests on a 
subset of this data against a program 30 that correctly deter- 
mines the errors, we found that the parameter values were 
consistent, and that the Igor errors exceeded the appropriate, 
one standard deviation (la) errors. In the following, our 
quoted parameter errors are those given by the program and 
should be considered the ler error. 

In analyzing the heat capacity data shown in Table II, we 
use the equations 

^ ^ + Et + A\t\ a +Z)|/| A| a (one phase), (6a) 

Cp~ + Et + A '|/| a + D'|/| Al a (two phase), 

(6b) 

where t = (T-T c )T c . In particular, we force equal exponent 
values above and below the critical point as predicted by 
theory. A smooth, continuous background specific heat near 
the critical point is expected, 3,6 which means the background 
amplitudes C p q and E are the same in the one- and two- 



Temperatirre CC) 

FIG. 2. The heat capacity C p of the fluid mixture triethylamine and water as 
given in Table II. Run 1 (O) and run 2 (X) are shown without their errors for 
clarity. A slight offset in the background heat capacity between the two runs 
is evident far from the critical point, which is at T t ~ 1 8.226±0.002. 


phase region. We simultaneously fit these equations to all the 
one- and two-phase data for each run. The exponent A, is 
fixed at its theoretical value of 0.5. This still leaves eight 
adjustable parameters: C pg , E, A, A \ D, D', T c , and a. The 
parameter E should reflect the noncritical fluid behavior and 
has a small negative value of -4 J/K for water in this tem- 
perature region which is effectively zero for the precision of 
our data. We do not find E to improve the fits (see Table III) 
and hold its value at zero for most fits of the heat capacity 
data. Allowing the critical exponent a to vary gives a value 
of 0.10±0.01, consistent with the theoretical value of 0.1 1, 
but substantially increases the errors in the amplitudes. Thus, 
most fits have a fixed. 

The two data runs reported in Table II and shown in Fig. 
2 have a small offset that is reflected in the value of C p . The 
critical temperatures, as determined by the fits shown in 
Table III, were slightly larger (about 2 mK) for run 2 than 
run 1 but not significantly so. The two runs give a good sense 
of our reproducibility. Fits 1-1 through 1-5 allow various 
parameters to vary or be fixed at the values indicated. Allow- 
ing a correction to scaling term (D,D‘ not zero) significantly 
improved the fit, but this form for our data does not allow a 



FIG 3. A plot of the enthalpy H values divided by the critical temperature 
T, for run 2 shows the small cusp at the critical point. 
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TABLE III. The heat capacity data from Table II are analyzed using Eq. (6) with A,— 0.5. A nonlinear, weighted least-squares fit is simultaneously done in 
the one and two phase region for each run. Units on C Pq , E , A, A’ , D, and D' , are J/(cm 3 K). The parameter errors are the one standard deviation uncertainty 
given by the fitting program. No error on a parameter means it was fixed. Values in brackets [ ] are calculated by imposing theoretical amplitude relations. 
The number before the dash in the “fit #’’ indicates the run. An * marks the best fit in our judgment. 


Fit # 

C po 

E 

A 

A' 

Tc 

D 

D f 

a 

r/N 

1-1 

1.13±0.03 

0 

1.78±0.02 

3.16±0.02 

18.2253 + 0.0006 

0 

0 

0.11 

2.06 

*1-2 

1.13±0.13 

0 

1 .78 ±0.05 

3.23±0.05 

18.2245 ±0.0006 

0.16±0.28 

- ].2±0.5 

0.11 

1.65 

1-3 

0.97 ±0.05 

-6.2± 1 .3 

l.86±0.02 

3.26±0.02 

18. 225 3 ±0.0006 

0 

0 

0.11 

1.72 

1-4 

0.4 ±0.7 

0 

2.3 ±0.5 

4.0±0.6 

18.225 ±0.001 

o 

1+ 

o 

Ui 

— 1.9±0.8 

0. 10±0.01 

1.61 

1-5 

1 . 14±0.26 

-4.8±2.4 

1.73±0. 18 

3. 10±0.23 

18.2251 ±0.0007 

0 

0 

0. 1 14±0.006 

1.74 

1-6 

1.37 ±0.1 3 

0 

1.69±0.05 

3.09±0.04 

18.2248 ± 0.0007 

-0.05±0.3 

[-0.9] 

0.11 

2.05 

1-7 

I.92±0.02 

0 

1.48±0.01 

2.90±0.01 

18.2238± 0.0005 

[-1.7] 

[-2.8] 

0.1 1 

2.53 

1-8 

1 .62 ±0.0 1 

0 

1.67±0.01 

[3.08] 

1 8.2265 ± 0.0007 

[-2.0] 

[-4.0] 

0.11 

7.97 

2-1 

0.79±0.03 

0 

1.94+0.02 

3.30+0.02 

)8.2278±0.0006 

0 

0 

0.11 

2.21 

*2-2 

0.96 ±0.1 3 

0 

I.87±0.05 

3.29±0.05 

18. 2269 ±0.0007 

-0.2+03 

— 1.510.5 

0.11 

1.94 

2-3 

1 . 1 5 ±0. 1 3 

0 

1.80±0.05 

3. 19±0.04 

18.2273 ±0.0007 

-0.8±0.2 

[-1.3] 

0.1 1 

2.13 

2-4 

1.69±0.02 

0 

1.58±0.01 

3.01 ±0.01 

18.2262+0.0007 

[—1.9] 

[-3.4] 

0.11 

2.41 

2-5 

1.46±0.01 

0 

1 .67 ±0.0 1 

[3.08] 

18. 2266 ±0.0007 

— 1.4±0. 1 

[-2.5] 

0.1 1 

2.19 


good determination of either D or D\ The best fit for run 1 
is fit 1-2 (see Table III) not only because the reduced chi 
square is lowest, but the resulting equation is the simplest 
that fits the data. Fits 1-6 through 1-8 impose amplitude 
relations 6,10,11 to reduce the number of variable parameters. 
Fit 1-6 allowed A, A \ and D to be free but calculated D ' 
assuming D + =D“. The theoretical value of R* was used 

°cr 

to calculate D and D' in fits 1-7 and 1-8 using the param- 
eters and A r . Fit 1-8 also imposes the theoretical 

constraint A + /A~=A/A f = 0.541, which allows one to fit the 
data with three parameters (C Pq ,A,T c ) but gives a poor fit to 
the data. It should also be noted that the definition of 

®cr 

allows one to choose the sign of D and we found negative 
values to give much better results. Bagnuls and Bervillier 13 
have argued that negative D values are acceptable theoreti- 
cally. None of these theoretical constraints designed to re- 
duce the number of variable parameters provided a good fit 
of run 1. 

The fits to the second run of the heat capacity data are 
also shown in Table III. Since run 2 is very similar to run 1, 
the parameter values are quite consistent. Corresponding fits 
for the two runs are the pairs: 2-1 and 1-1, 2-2 and 1-2 (the 
best fits), 2-3 and 1-6, and 2-4 and 1-7. Fit 2-5 shows that a 
fairly good fit can result when A/ A ' is forced to be 0.541 and 
D + 1 D~ is forced to be 1.0. This fit can be used to calculate 
a value of Rg = —0.63, significantly less than the 
predicted 11 value of -0.71. 

Fits can also be done on the enthalpy, which is calcu- 
lated by integrating the net power into the fluids over time as 
in Eq. (4). Using Eq. (5) with the expression for the heat 
capacity given in Eq. (6), we can express the enthalpy H as 


H H c E A 

T c ~ T c + Cp o l+ 2 ' + ( 1 - a) ^ 


D 

( I + A, — a) 


M 


I + A | * a 


(one phase), 


(7a) 


H 

T~ 

1 c 


H c E A’ 

T c +Cp o t+ 2 ' + (l-a) 


M' 


D' 

+ (1 +A,-a) 


, + A| a (two phase), 


(7b) 


where t is the reduced temperature, H c is the enthalpy at the 
critical point, and the other symbols are the same as in Eq. 
(6). The enthalpy has a cusp (1-a) anomaly at the critical 
point which is difficult to resolve when compared to the 
background (noncritical) quadratic dependence. Moreover, 
by integrating the raw data we sum any residual systematic 
errors in the experiment which appear as part of the “back- 
ground” enthalpy. Since we take into account all known 
systematic errors in calculating the enthalpy, each data point 
is assigned the same error of lXKT 5 J/(cm 3 K) in HfT l . 
With the large number of data points (10 274 for run I and 
1 1 643 points for run 2), most of which are far from critical, 
the fits are inherently weighted to fit the region away from 
critical. Nevertheless, the resulting parameter values are 
quite consistent with those determined by the heat capacity 
fits (see Table IV). 

In order to better determine the critical behavior, the 
range in temperature was restricted (| r | <2X 1 0 -3 ) for some 
of the fits to the enthalpy for each run. By using enthalpy 
values close to the critical point, the correction to scaling 
terms D and D ' should be zero while the cusp anomaly with 
amplitudes A and A ’ should be significant. The exponent a 
was free in some of the fits and its value was determined to 
be 0.107±0.006, which is quite consistent with theoretical 
predictions 9 of 0.1 10±0.003. We fix the value of a for most 
of the fits so that the amplitudes, which are sensitive to the 
exponent value, can be compared to other experimental and 
theoretical determinations where a was set at 0. 1 1 . The re- 
duced range fits 1-D and 2-D in Table IV show that the 
amplitudes A and A ' are unchanged from the full range fits 
1 -B and 2 -A in which correction to scaling terms are in- 
cluded. Moreover, the reduced range fits I -E and 2-F show 
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TABLE IV. The enthalpy H , as calculated by Eq. (4) and in units of J/cm 3 is fitted hv Fn (1\ , . 

critical temperature was fixed at 18.225 for run 1 and 18.227 for run 2 ’ y Eq. ( ). The parameters and notations are defined in Table III. The 



1-A 

*l-B 

1-C 

1- D 
1 -£ 

*2-A 

2- B 
2-C 
2-D 
2-E 
2-F 


full 
full 
full 

M<2X10 -3 

j/| <2 X 10 -3 

full 

full 

full 

|/|<2XI0~ 3 
|f|<2XI0~ 3 
|rj <2x 10 -3 


5.544 ±0.001 
5.538±0.001 
5.541 ±0.001 
5.541 ±0.001 
5.538±0.001 
6.215±0.001 
6.215±0.001 
6.216±0.001 
6.215±0.001 
6.215±0.001 
6.214±0.001 


~p o 

U9±0.01 
1.53±0.01 
1.45±0.12 
1.42 ±0.01 
1.45±0.01 
1.09 ±0.01 
0.97 ±0.03 
1. 12 ±0.04 
l.I6±0.01 
1.31 ±0.03 
1 - 1 5 ±0.0 1 


: 1 


0 
0 

-19 
0 

-54±1 

0 

0 

- 11.8 + 0.1 
0 
0 

-26.7±0.2 


1.787±0.002 
1. 634 ±0.006 
1.84 ±0.09 
1.674 ±0.006 
1.634 ±0.006 
1.838 ±0.002 
1.94 ±0.02 
1.94 ±0.03 
1.797 ±0.006 
1.69 ±0.02 
1.791 ±0.006 


3.179±0.002 

3.105±0.006 

3.39 ±0.12 
3.093 ±0.006 
3.105±0.006 
3.267 ±0.006 

3.40 ±0.03 

3.41 ±0.03 
3.203 ±0.006 
3.06 ±0.02 
3.215±0.006 


0.11 
0.11 

0.101 ± 0.002 
0.11 
0.11 
0.11 

0. 107±0.001 
0.105±0.001 
0.31 

0.1 13±0.002 
0.11 


0 

-0.55 ±0.03 
~2.4±0.2 
0 
0 

— 0.46±0.03 
-0.4±0.1 
— 1.9±0.I 
0 
0 
0 


-2.33±0.07 
— 3.0±0. 1 
0 
0 

-1.81 ±0.07 
- 2 . 0 ± 0. 1 
-2.4±0.1 
0 
0 
0 


4.4 

2.04 
1.67 
4.15 

2.95 
2.31 
2.27 

1.95 
2.00 
1.97 
1.69 


that including E =£ 0 improves the fit but does not change the 
values of the amplitudes A and A The same is not the case 
for the full range fits 1-C and 2-C where allowing E^0 
causes the amplitudes A and A ' to increase. Therefore, we 
are more confident in the consistent values of A and A' de- 
termined when E=0. The enthalpy data require the use of 
extended scaling for the entire range (|r|<10~ 2 ) but the val- 
ues depend heavily on the whether E is zero or free (compare 
fits 1 -B to 1-C and 2 -A to 2-C). Because of the consistency 
of A and A' in the restricted and full range fits, we favor the 
values of D and D' as determined in fits 1 -B and 2-A, which 
is why we label these as our best fits over the entire range. 

The heat capacity fits given in Table III and the enthalpy 
fits in Table IV provide a consistent set of amplitude values 


when E 0 and a— 0.1 1, which can be converted to the am- 
plitudes of Eq. (1) to be: £ c . = -2.7±0.2, A + =0.201 ±0.006, 
A =035S±0.005 J/(cm 3 K), Z> + = -0.29 ±0.04, and 

D =-0.65±0.10, where all the errors are one standard de- 
viation estimates. This value of A + is somewhat larger, while 
the value of A is somewhat smaller, than Thoen et aO 5 
found in their heat capacity experiment, yet both are within 
experimental error of those earlier values (see Table I). How- 
ever, A is inconsistent with all the values determined using 
either density or refractive index measurements, whose val- 
ues tend to be inconsistent with each other. Our value of the 
universal ratio A + /A “ =0.56±0.02 was quite consistent with 
experimental values determined in several binary-liquid sys- 
tems recently (see Table V) and slightly larger than theoreti- 


I- L ^p V a Amphiu ** m[os dimmed from this experiment. Various techniques have been used to determine 
A in TEA- water (values in Table I) and determine X using assorted values of our determination uses the 
r^strecem y measured value. Three other recently measured binary liquid systems (triethylamine-heavy water 
[TEA-D 2 0], 3-methy Ipentane 4- nitroethane [3MP-NE], and 2,6-lutidine+ water [2,6 L-H 2 0]) are provided for 

srm.r™* p " ai “ d 


System/technique 

A*IA~ 

D + /D~ 

X 

R l, 


tea-h 2 o 

p 

n 

C r 

tea~d 2 o 
3MP-NE 
2,6 L-H 2 0 

0.36±0.03 a 
0.71 ±0.20 b 

0.56±0.02 d 

0.57±0.01 f 

0.56±0.02 g 

0.57±0.25 h 

0.5±0.1 d 

— 0.5 C 
I.0±0.4 g 

0.01 7 ± 0.00 l a 
0.018±0.002 a 
0.017±0.006 b 
0.024 ±0.007 c 

0.019± 0.003 d 

0.019±0.004 e 
0.020 ±0.002 g 

-0.5±0.1 d 

— 0.48±0.01 f 

-u±0.1 d 

-0.96±0.03 f 

Liquid-vapor 

0.52 ±0.02' 

1.4* 




Theory 

series 

RG 

0.523±0.009 k 
0.541 ±0.014 m 

0.96±0.25 m 

0.0188±0.0002 k 
0.0 197 ± 0.0002" 

-0.708 15 p 

- 1.334±0.044 m 

Reference 27. 







Reference 15. 

Reference 21. 
d Our results in bold. 

'Reference 17. 

Calculated from parameter values 
in Table III of Ref. 17. 
Reference 18. 


Reference 19. 
Reference 2. 
k Reference 9. 
"Reference 6. 
"Reference 14. 
Reference 1 1 . 
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cal predictions using either renormalization group 6 or series 9 
calculations. The amplitude ratio D^ID " =0.5 ±0.1 is no 
longer thought to be universal 13 but is expected to be similar 
in similar systems. The value we find is the same magnitude, 
but opposite sign, of that found in the triethylamine-heavy 
water system, 17 and about one half the value found in 
3-methylpentane + nitroethane ls and predicted by theory. 6 Fit 
2-5 indicated that our data are not inconsistent with 
D^ID =1, but these difficult to determine, correction to 
scaling amplitudes are best determined in our enthalpy tits, 
which favor a lower ratio. 

Two other universal amplitude ratios can be determined 
from our experiment. The two-scale-factor universality ratio 
X = A } can be calculated from our value o (A f and the 
most recently published 21 value for §, = 0. 1 10 ±0.005 nm, 
which gives X=0.019±0.003. This value agrees with other 
recent experimental determinations and with theoretical pre- 
dictions (see Table V). 

It is not certain whether another amplitude ratio R # 
given by Eq. (2) is universal; 7 S l ~ 13 however, theoretical pre- 
dictions have been made 6 11 which have not previously been 
tested in a binary fluid mixture. Even though we could not 
determine D with much precision, it enters into this ampli- 
tude ratio raised to a small power and thus is not so crucial. 
We find R f \ = — 0.5 ± 0. 1 and R H ~ — I . I ± 0. 1 , w hich are 
both smaller than predicted (see Table V). Using the fit pa- 
rameters Bloemen el al. ]l determined for the system triethy- 
lamine and heavy waiter, we can calculate values for R f{ 
which are consistent with ours (see Table V). Belyakov and 
Kiselev 8 feel that it would be difficult to distinguish the sin- 
gular specific heat B { from the regular part and thus making 
an experimental check of the universality of /?„ “probably 
impossible.” The system studied here has the advantage of a 
constant specific heat in the one-phase region far from the 
critical temperature [£ /j( ,=4.07 J/(cm 3 K) in run 1 and 4.01 
J/(cm 3 K) in run 2], which together with the fitted C/> allows 
a good approximation to B ( . It is thus possible to begin the 
process of experimental determination of a possibly univer- 
sal value for R ^ . 

CONCLUSION 

The heat capacity of the binary liquid mixture 
triethylamine-water has been measured near its lower criti- 
cal consolute point using a scanning, adiabatic calorimeter. 
The critical exponent a was determined to be 0.1 07 ±0.006, 
consistent with theoretical predictions. While our values of 
A f =0.201 ±0.006 and A =0.358 ±0.005 J/(cm 3 K) for 
<*-0.11 were consistent with a previous heat capacity 
measurement, 15 the value of A h w'as inconsistent w ith values 
determined by density or refractive index 
measurements- 6 *" 0- " 3 in this triethylamine and water system. 
Jacobs and Greer 31 have reviewed the values of A " as ob- 
tained from density measurements and find them not consis- 
tent with those obtained from C fJ measurements; they sug- 
gest that the discrepancy is due to ignoring the behavior of 

^ T’ ^ 

c v near l ( . - 


Various amplitude ratios were calculated from the best 
fits to our data with a=0.11: A +IA ~ =0.56±0.02, 
D~ ID =0.5±0. 1, X = 0.0 19 ±0.003, R^ = -0.5 ±0.1, and 
Rfi — —1.1 ± 0. 1 . While A '/A ~ was consistent with other 
recent experimental determinations in binary liquid mixtures, 
it was slightly larger than both theoretical predictions and 
recent experimental determinations in the simple liquid-gas 
systems C0 2 and SF 6 . 19 The correction to scaling amplitude 
ratio DID was half of that predicted, but experimental 
values are still quite scattered. As a result of more precise 
theoretical calculations and experimental determinations, the 
two-scale-factor universality ratio X now' is consistent among 
experiments and theories. We could not reduce the number 
of variable parameters in our fits by using R ^ , but that is 
probably because its universal value has not been deter- 
mined. Our determination of R^ is smaller than predicted 
and is the first such determination in a binary fluid mixture. 
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